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Abstract 

Recently, two measures of electroweak finetuning (EWFT) have been introduced for SUSY 

'~~^ models: Aew compares the Z mass to each separate weak scale contribution to m^z 

O, while Ahs compares the Z mass to high scale input parameters and their consequent 

' renormalizaton group evolution (1/A is the % of fine tuning). While the paradigm 

Q mSUGRA/CMSSM model has been shown to be highly finetuned under both parameters 

r^i {Aew ~ 10^ and Ahs ~ 10^), the two-parameter non-universal Higgs model (NUHM2) 

in the context of radiatively-driven natural SUSY (RNS) enjoys Aew as low as 10, while 

►^ Ahs remains ~ 10^. We investigate finetuning in the 19- free-parameter SUGRA model 

C\| (SUGRA19). We find that with 19 free parameters, the lowest Aew points are comparable 

^^ to what can be achieved in NUHM2 with just 6 free parameters. However, in SUGRA19, 

Ahs can now also reach as low as 5 — 10. The conditions which lead to low Ahs include 



iTT-Hu ~ ^^z at the high scale, with non-universal gaugino masses Mi^2 S> M^ also at 

(O niGUT- The low Ahs rnodels are severely constrained by 6 — )• ,57 branching fraction. In 

CO both cases of low Aew and Ahs, the superpotential n parameter should be ~ 100 — 300 

. . GeV. While SUSY models with low EWFT may or may not be discoverable at LHC, the 

.!^ predicted light higgsinos must show up at a linear e+e^ collider with ^/s ~ 2|//|. 
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1 Introduction 

Supersyinmetric models of particle physics are renown for providing an elegant solution to the 
daunting gauge hierarchy problem. They also receive indirect experimental support from 1. 
the measured strengths of weak scale gauge couplings, which allow for unification at a scale 
itlgut — 2 X 10^^ GeV within the Minimal Supersymmetric Standard Model (MSSM) and 2. 
from the measured value of the top quark mass, which is sufficiently high as to radiatively 
drive electroweak symmetry breaking (REWSB)[lJ. Along with these plaudits, 3. the recent 
discovery by Atlas [2] and CMS [3] of a Higgs-like boson with mass mh — 125 GeV confirms 
predictions from models of weak scale supersymmetry|l] where (in the context of the MSSM) 
a value m^ ~ 114 — 135 GeV was required p]. The emergent picture is that the MSSM (or 
possible extensions) may provide a solid description of nature not only at the weak scale, but 
perhaps all the way up to energy scales associated with grand unification p]- 

Such an audacious extrapolation has suffered a string of serious set-backs: so far, no signs of 
supersymmetric matter have emerged from LEP, LEP2, Tevatron or, more recently, LHC data. 
Recent analyses from Atlas and CMS in the context of the minimal supergravity (mSUGRA 
or CMSSM) model[7] require mg ~ 1.4 TeV for mq ~ nig and mg ~ 1 TeV for mg <^ mq. 
Naively, these results exacerbate the so-called Little Hierarchy Problem (LHP) [8] : why is there 
an apparent discrepancy between the weak scale (typified by mz — 91.2 GeV) and the SUSY 
scale, where m{sparticle) ~ 1 TeV? The growing scale mismatch has led some physicists to 
call into question whether or not weak scale SUSY really exists, or at least to concede that it 
suffers unattractive electroweak finetunings (EWFT) [H] . 

Traditionally, EWFT has been quantified using the Barbieri-Giudice measure [TI H [TT | [T ^ [T3 | 
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where Oj represents various fundamental parameters of the theory, usually taken to be some 
set of soft SUSY breaking parameters defined at some high energy scale Khs below which the 
theory in question is posited to be the correct effective field theory description of nature. 1/A is 
the % of fine tuning. The value of A^g then answers the question: how stable is the fractional 
Z-boson mass against fractional variation of high scale model parameters? Depending on which 
parameters are included in the set Oj, very different answers emerge [21]. In addition, theories 
which are defined at very different values of Is^hsi but which nonetheless lead to exactly the 
same weak scale sparticle mass spectra, give rise to very different values of A^^. 

1.1 AijvF and A^^^ 

Recently, two different measures of EWFT- ^ew and ^hs^ have been proposed which answer 
a different but related question: how is it possible that mz has a value of just 91.2 GeV while 
gluino and squark masses exist at TeV or even far beyond values? The answer should be: 
those independent contributions which enter the scalar potential and conspire to build up the 
Z-boson mass should all be comparable to niz- 



Minimization of the scalar potential in the MSSM[1] leads to the well-known relation that 

m| _ mj,^ + S^ - (m^^ + E-) tan^ /3 ^ ..^ 

2 " tan2/3-l ^ ' ^ ' 

where fn'jj^ and m|^ are soft SUSY breaking (not physical) Higgs mass terms, /i is the super- 
potential Higgsino mass term, tan/3 = Vu/v^ is the ratio of Higgs field vevs and S^ and SJ^ 
include a variety of independent radiative corrections |25]. 



1.1.1 A 



EW 



Noting that all entries in Eq. [2] are defined at the weak scale, the electroweak fine-tuning 
parameter 

Aew = maxi \Ci\ /(m|/2) , (3) 

may be constructed, where Ch^ = mj^^/ (ta.n^ (3 — 1), Ch^ = — m|^^ tan^/3/(tan^ /3 — 1) and 
C^ = -/i^. Also, Csu(fc) = S^(A;)/(m|/2) and Cs^(fc) = ^dik)/{m%/2), where k labels the 
various loop contributions included in Eq. |2} A low value of Aew means less fine-tuning, 
e.g. Aew = 20 corresponds to A^J^ = 5% finetuning amongst terms contributing to m|/2. 
Since Ch^ and (^^^(fc) terms are suppressed by tan^/3 — 1, for even moderate tan/3 values the 
expression Eq. [2] reduces approximately to 



rrir. 
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In order to achieve low Aew, it is necessary that — m|^ , —fi"^ and each contribution to — S^ all 
be nearby to Tn\/2 to within a factor of a few. 

A scan over mSUGRA/CMSSM parameter space, requiring that LHC sparticle mass con- 
straints and iTLh = 125 ± 2 GeV be obeyed, finds a minimal value of Aew ~ 10^, with more 
common values being Aew ~ 10^ — 10^. Thus, one may conclude that the Z mass is rather 
highly finetuned in this paradigm model. In the case of mSUGRA, the value C^ becomes low 
only in the hyperbolic branch/focus point [121 [IB] (HB/FP) region. In this region, however, tuq 
and consequently rn^^^ are very large, so that S^(ti^2) are each large, and the model remains 
finetuned. 

Alternatively, if one moves to the two-parameter non-universal Higgs model (NUHM2) |26] , 
with free parameters 

mo, mi/2, Aq, tan/3, /i, m^ (5) 

then 

1. fi can be chosen in the 100 — 300 GeV range since it is now a free input parameter, 

2. a value of m?u^{mGUT) ~ (1-3 — 2.5)?7io may be chosen so that Tn\j^ is driven only slightly 
negative at the weak scale, leading to m\j {weak) ~ — ?7i|/2, and 

3. with large stop mixing from Aq ~ ±1.6?7io, the top-squark radiative corrections are soft- 
ened while rrih is raised to the ~ 125 GeV level |25j. 



In the NUHM2 model, /S.ew as low as 5 — 10 can be generated. For such cases, the Little 
Hierarchy Problem seems to disappear. The low ^ew rnodels are typified by the presence of 
light higgsinos m~ , m^ ~ 100 — 300 GeV which should be accessible to a linear e"'"e~ collider 

operating with a/s ~ 2|yu|. Also, nig ~ 1 — 5 TeV while mi_^ ~ 1 — 2 TeV and m^^ ~ 2 — 4 TeV. 
The measure ^ew listed above is created from weak scale MSSM parameters and so contains 
no information about any possible high scale origin, even though low values of Aew rnay be 
required of high scale models: in this sense, low /S.ew captures a minimal EWFT required of 
even high scale SUSY models. 

1.1.2 Ahs 

To include explicit dependence on the high scale A at which the SUSY theory may be defined, 
we may write the weak scale parameters m%^ ^ and /i^ in Eq. (2| as 

"^L,. = ^L,. (^) + '^^L,. 5 /"' = /"' (A) + 5l^^ , (6) 

where m\j (A) and ^Ji^{h) are the corresponding parameters renormalized at the high scale 
A. It is the (5m|^ terms that will contain the log A dependence emphasized in constructs of 
natural SUSY models pTf [251 I2U] . In this way, we write 

ml (m?. (A) + (5m2 + S^) - (m2 (A) + (5m2 + S") tan2 /3 ^„^^ „ 

z tan p — 1 

In the same spirit used to construct Aew, we can now define a fine-tuning measure that encodes 
the information about the high scale origin of the parameters by requiring that each of the terms 
on the right-hand-side of Eq. Q (normalized again to m|/2) be smaller than a value Ahs. The 
high scale fine-tuning measure Ahs is thus defined to be 

Ans = maxi\Bi\/{m%/2) , (8) 

with Bhj^ = m|^^(A)/(tan^ /? — 1) etc., defined analogously to the set Cj. 

As discussed above, in models such as mSUGRA whose domain of validity extends to very 
high scales, because of the large logarithms one would expect that (barring seemingly accidental 
cancellations) the BsHu contributions to Ahs would be much larger than any contributions to 
Aew because the term m'jj^ evolves from large ttiq through zero to negative values in order to 
radiatively break electroweak symmetry. Thus, Ahs is numerically very similar to the EWFT 
measure advocated by Kitano-Nomura^T] where A^n = <^"^ii-„/("^^/2) 

Scans of the mSUGRA/CMSSM model in Ref. [3D] found Ahs ~ 10^. In Ref. [2S], scans 
over NUHM2 model similarly found Ahs ~ 10^. Thus, both the mSUGRA and NUHM2 models 
would qualify as highly EW finetuned under Ahs- 

1.2 Goals of this paper 

In this paper, we would like to maintain the SUSY grand desert scenario where the MSSM is 
postulated as the correct effective theory below Q ~ mcuT- However, we would like to expand 



our set of input parameters, in this case, to a maximal set of 19, which maintains the scenario of 
minimal flavor and minimal CP-violation. The resulting model, dubbed here as SUGRAl9[31j. 
has the same parameter freedom as the more popular pMSSM model[32]. However, unlike 
pMSSM defined at the weak scale, SUGRA19 maintains the successes of renormalization group 
evolution, and its consequent gauge coupling unification and radiative electroweak symmetry 
breaking due to the large value of rrit. 

In this paper, we have several goals. The first is to check, under models with maximal 
parameter freedom, whether even lower values of Aew ^ 10 can be found, or whether NUHM2 
already achieves the minimal EWFT values. We will find that Aew ~ 5 — 10 is about as low 
as can be achieved while maintaining accord with phenomenological constraints, and that the 
resulting models tend to look phenomenologically rather similar to RNS models as derived from 
NUHM2 parameter space. 

Our second goal is to check whether values of Ahs ^ 10^ can be found. In the case of 
SUGRA19, we will find that Ahs as low as 5 — 10 can also be found, but only for special choices 
of non-universal SUGRA parameters. The low Ahs models appear tightly constrained if accord 
with BF{b — )■ s'j) measurements is imposed. Both the low Aew and the low Ahs niodels are 
characterized by light higgsinos with mass ~ 100 — 300 GeV, which should be accessible to 
linear collider searches, although perhaps not accessible to LHC searches. We present some 
interesting benchmark (BM) points for both low Aew and low Ahs- 

Before proceeding, we mention our results in relation to several previous works on reduced 
EWFT in models with non-universal soft SUSY breaking terms. After several initial studies 
using the BG measure in universal models [TUl IH [121 [131 [U], Kane and King[17j showed that 
non-universal gaugino masses with Mi, M2 > M3 lead to reduced EWFT. These were followed 
by similar studies by Lebedev et a/. [19], Abe et a/. [20], Martin [21] (the latter of which was used 
to motivate a compressed SUSY mass spectrum) and Antusch et a/. [22]. Recently, Gogoladze 
et al.[33\ studied the measures Aew and Ahs within the context of GUT-motivated gaugino 
mass non-universality using effectively a 5-parameter model with gaugino masses related by 
Ml = IM3 + |M2. In their study, they were already able to reduce the maximum Ahs values 
down to the ~ 30 level, which is already quite close to what is achieved here using 19 free 
parameters. 

2 Scan over 19 parameter SUGRA model 

To calculate superparticle mass spectra in the SUGRA19 model, we employ the Isajet 7.83 [34j 
SUSY spectrum generator Isasugra[35]. Isasugra begins the calculation of the sparticle mass 
spectrum with input DR gauge couphngs and fb, fr Yukawa couplings at the scale Q = Mz 
{ft running begins at Q = rrit) and evolves the 6 couplings up in energy to scale Q = Mqut 
(defined as the value Q where gi = (72) using two-loop RGEs. We do not enforce the exact 
unification condition gs = gi = g2 at Mqut, since a few percent deviation from unification 
can be attributed to unknown GUT-scale threshold corrections [36] ■ Next, we impose the 
SSB boundary conditions at Q = Mqut and evolve the set of 26 coupled two- loop MSSM 
RGEs ^7\ [3S] back down in scale to Q = Mz- Full two-loop MSSM RGEs are used for soft 
term evolution, and the gauge and Yukawa coupling evolution includes threshold effects in 



the one-loop beta-functions, so the gauge and Yukawa couphngs transition smoothly from the 
MSSM to SM effective theories as different mass thresholds are passed. In Isasugra, the values 
of SSB terms which mix are frozen out at the scale Q = msusY = y/'T^^jW^^ while non-mixing 
SSB terms are frozen out at their own mass scale [35]. The scalar potential is minimized 
using the RG-improved one-loop MSSM effective potential evaluated at an optimized scale 
Q = iTT'SUSY to account for leading two-loop effects [39j. Once the tree- level sparticle mass 
spectrum is obtained, one-loop radiative corrections are calculated for all sparticle and Higgs 
boson masses, including complete one-loop weak scale threshold corrections for the top, bottom 
and tau masses at scale Q = nisusY [10] • Since Yukawa couplings are modified by the threshold 
corrections, the solution must be obtained iteratively, with successive up-down running until a 
convergence at the required level is found. 

We search for models with low Aew and low A^s by first performing a broad-based random 
scan over the following SUGRA19 parameter ranges: 

• Gaugino masses: Mi, M2, M3 : - 3.5 TeV 

• First/second generation scalar masses: mq-^, mu^, mjj^, rriLi, ttie^: — 3.5 TeV, 

• Third generation scalar masses: mq^, f^Us, '"^Ds, f^Ls, 'rriE^'- — 3.5 TeV, 

• Higgs soft masses: rnn^, ^^h^ '■ ^ 3.5 TeV, 

• trilinear soft terms: At, A^, A,-:— 3.5 TeV — )■ 3.5 TeV, 

• ratio of weak scale Higgs vevs tan/3 : 2 — 60. 

We adopt a common mass for first and second generation scalars so as to avoid the most 
stringent SUSY FCNC constraints [41j. 
We require of our solutions that: 

• electroweak symmetry be radiatively broken (REWSB), 

• the neutralino Zi is the lightest MSSM particle, 

• the light chargino mass obeys the model independent LEP2 limit, m^ > 103.5 GeV|42j 
and 

• 123 <mh< 128 GeV. 

We do not impose any LHC sparticle search limits since our general scan can produce com- 
pressed spectra which in many cases can easily elude LHC gluino and squark searches. Points 
which satisfy the above constraints are plotted as blue circles in the following scatter plots. 

We will also calculate the values of BF{b -^ 57) |13l IB] and BF{Bs -)■ /i"^/x")[lS] for each 
point generated. The measured value of BF{b — )■ 57) is found to be (3.55 ± 0.26) x 10~^ [3S]. 
For comparison, the SM prediction [17] is BF^^'ib -)■ S7) = (3.15 ±0.23) x 10"^. Also, recently 
the LHCb collaboration has found an excess over the background for the decay Bg — )■ /i"*"/^" [H] . 
They find a branching fraction of BF{Bs — >■ ft~^fJ'~) = S.2t\'^ x 10~^ which is in accord with the 
SM prediction of (3.2 ± 0.2) x 10^^. Points with BF{b — )■ 57) within 3cr of its measured value 
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Figure 1: Plot of Ahs vs. Aew from a broad (dark/light blue) and focused (red/orange) 
scan over SUGRA19 model parameter space. The orange and light blue points satisfy -B-decay 
constraints while the dark blue and red points do not. 



BF{b -^ S7) = (2.5 - 4.5) x 10"^ and points with BF{Bs -^ /i+/i") = (2 - 4.7) x lO^^ will be 
labeled as light blue, showing that these points are also in accord with S-physics constraints. 

Our first set of results are shown in Fig. [1} The broad scan points are shown in blue. We 
see that the bulk of generated points yield Aew and Ahs ~ 10^, so would qualify as highly EW 
finetuned in generating mz = 91.2 GeV. The points with the lowest Aew values come in with 
Aew ~ 10, which is similar to that which can be achieved in the more restrictive NUHM2, but 
which is much better than what can be achieved in mSUGRA. 

The lowest Aew point has Aew = 7.9, while the corresponding A^s = 190. The SUGRA19 
parameters associated with this point are listed in Table [T] in the column labeled as EWl. The 
point has the required low /i ~ 180 GeV and m?H^{mweak) ~ —(171 GeV)^. In addition, the 
large top-squark mixing At/mQ{?)) ~ —2.1 softens the top squark radiative corrections S^(ti^2) 
whilst raising rrih up to 123.5 GeV. 

The corresponding sparticle mass spectra are listed in Table 12} The gluinos and squarks 
are ~ 2 — 3 TeV: well beyond current LHC reach. The W2 and Z12 are dominantly higgsino- 
like with a mass gap m^ — m^ ~ 3 GeV. Thus, even though the higgsinos can be produced 
with large cross sections at LHC, the very soft visible energy release from their decays makes 
them difficult to detect |19]. The light higgsinos should be straightforward to detect at a linear 
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parameter 


EWl 


HSl 


HS2 


MiimGur) 


2822.1 


3266.2 


3416.4 


M2{mGUT) 


3385.3 


2917.8 


3091.3 


M^{mGUT) 


884.9 


1095.7 


1085.8 


mQ{l) 


2484.7 


1192.6 


978.5 


mu{l) 


2506.2 


2468.3 


2440.6 


rriDiX) 


2342.1 


1508.9 


1404.2 


mi(l) 


1820.4 


623.8 


754.8 


mE{l) 


1731.2 


936.1 


915.8 


^q(3) 


698.3 


6.6 


371.3 


mc/(3) 


1552.8 


233.9 


23.2 


mz)(3) 


1498.5 


2946.0 


3052.4 


mL(3) 


3339.3 


341.1 


451.3 


m£;(3) 


2114.9 


1268.7 


1247.5 


"^H„ 


871.3 


314.0 


125.4 


^H^ 


2205.3 


3160.4 


2964.9 


At 


-1509.6 


-1024.4 


-801.3 


A, 


2301.7 


3121.6 


3294.3 


Ar 


3307.3 


1932.0 


1754.5 


tan/3 


27.0 


51.1 


29.0 


/^ 


181.4 


242.8 


98.0 


Aw 


7.9 


17.9 


5.2 


Ah5 


190.0 


32.0 


6.4 



Table 1: Input parameters (GUT scale) in GeV for one low /S^ew point and two low A//5 points. 
We take rrit = 173.2 GeV. 

e~^e~ collider with ^/s ~ 400 GeV. The lightest top squark ti has mass less than 1 TeV: 
typically below values generated from radiative natural SUSY models[25]. This leads to a 
somewhat anomalous branching fraction BF{b — )■ S7) ~ 2.5 x 10"'^, below the measured value 
of (3.55 ± 0.26) X 10"^ gS]. 

A perhaps surprising result from Fig. IT] is that A^s values far below the NUHM2/mSUGRA 
minimal value of 10^ can now be found. In fact, the lowest A^s point from the broad scan has 
a value of 32, or 3.1% EWFT, even including the effect of high scale logarithms. The parameter 
values for this point, labeled as HSl, are also hsted in Table [l] There are several features of the 
input parameters which lead to low Ahs- First, the GUT scale value of m|^^ = (314 GeV)^, so 
our high scale starting point for m/f^ is not too far from ruz- Second, the GUT scale gaugino 
masses Mi and M2 are ~ 3M3 ~ 3 TeV. The RG running of mfj is governed by 

'^'"^"- ^ ^-%M-39lM'2 + ^9!S + 3f!Xt) (9) 



dt 167r2 V 5''' ' ""^ ^ 10^ 

where t = log((5^/yU^), 5* = rajj^—mfj^ + Tr rng — m| — 2m^ + m|, + m|; and Xf = mQ(3) + 
mfj{3) + rn'jj^ + Af. At Q = ttigut, the large gaugino masses provide a large negative slope 
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HSl 


HS2 


m~g 




2042.9 


2436.7 


2428.8 


rriu^ 




3650.7 


2991.9 


2968.5 


mu^ 




2980.5 


3214.8 


3191.6 


^en 




2196.3 


1763.6 


1786.1 


mi 




879.5 


1033.2 


892.4 


^2 




2305.1 


1958.3 


2394.9 


'T^k 




2121.8 


1961.4 


2418.0 


T^h 




2327.7 


2916.1 


3495.8 


rrifi 




2219.6 


1049.5 


1748.3 


mf2 




3865.8 


1467.5 


1911.3 


ruy^ 




3884.8 


1464.9 


1911.4 


m^. 




2802.2 


2393.0 


2538.3 


"^m 




192.1 


255.5 


104.1 


^z. 




2810.2 


2386.8 


2530.3 


"^Z3 




1261.2 


1448.0 


1513.5 


"^z. 




187.8 


251.2 


102.4 


^z. 




184.7 


247.9 


99.3 


rriA 




2759.7 


2242.6 


3176.4 


rrih 




123.5 


123.6 


123.1 


nfK' 




0.007 


0.013 


0.003 


BF{h 


-^ S7) X 10^ 


2.5 


1.8 


2.6 


bf\b 


s -> /i+yU") X 10*^ 


3.9 


4.5 


3.8 


cy'\Z, 


P) (Pb) 


2.9 X 10'^° 


3.7 X 10-1° 


2.5 X 10-1° 



Table 2: Sparticle masses in GeV and observables for one low IS.ew and two low l\us points as in 
Table^ The measured values of the branching fractions are BF{b — )■ s'j) = (3.55±0.26) x 10"'^ 
and BF{Bs -> yU+/i-) = 3.2+};^ x 10"^ 
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Figure 2: Plot of a) slope dm'jj^/dt vs. Q from model HSl with Ahs = 32. 



(green curve of Fig. pi) for mj^^ , causing its value to increase while running towards lower mass 
scales. As the parameters evolve, Xt increases due to the increasing squark soft terms so that 
the Yukawa coupling term grows (red curve from Fig. pi) and ultimately dominates; then mjj^ 
is driven towards negative values, so that electroweak symmetry is finally broken. The total 
slope (black curve) passes through zero around Q ~ 10^° GeV, indicating large cancellations 
in the RG running of Tn'jj^. Ultimately, the value of Tn?jj^{mweak) ~ —(185 GeV)^ so that both 
the starting and ending points of fn^j^ remain not too far from r7i|, and hence 8m\j^ is not too 
far from m|, thus fulfilling the most important condition required by low /S.hs- 

The RG running of gaugino masses and selected soft scalar masses for HSl are shown in 
Fig. p] In frame a), we see that indeed Mi and M2 start at ~ 3 TeV values and decrease, whilst 
M3 starts small at Q = uigut and sharply increases. The gaugino mass boundary conditions 
then influence the running of the soft scalar masses in frame h). Most important is the running 
of m?Hu^ which starts near m| at mcuT, Tuns up to about the TeV scale at Q ~ 10^° GeV, and 
then is pushed to small negative values by Q ~ i^weak- Also, m[/(3) and mqi?)) start small, 
which aides the high Q gaugino dominance in the running of rn^j^. By Q ~ rriy^eak, these third 
generation squark soft terms have been pushed to the TeV scale. Thus, top squarks are not so 
heavy and the radiative corrections E"(ti^2) are under control. 



3 Results from narrow scan 



To hone in on SUGRA19 solutions with low Ahs, we will impose a narrow, dedicated scan 
about our lowest Ahs solution: 



Ml : 3 - 3.5 TeV, M2 : 2.7 - 3.2 TeV, M3 : 0.8 - 1.3 TeV 

mQ(l,2) : 0.9 - 1.4 TeV, mc/(l,2) : 2.2 - 2.7 TeV, mD(l,2) 
mL(l,2): 0.4-0.9 TeV, mi5(l, 2) : 0.7 -1.2 TeV, 



1.25 - 1.75 TeV, 
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Figure 3: Plot of a) running gaugino masses and b) running scalar masses vs. Q from model 
HSl with Ahs = 32. 
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• mQ(3) : - 0.5 TeV, mt/(3) : - 0.5 TeV, mz?(3) : 2.7 - 3.2 TeV, mi(3) : 0.1-0.5 
TeV, niEiS): 1-2 TeV, 

• niH.^ : 0.05 - 0.55 TeV, nin^ : 2.9 - 3.4 TeV, 

• At : -1.3 -^ -0.8 TeV, A^ : 2.9 - 3.4 TeV, Ar : 1.7 - 2.2 TeV, 
with tan /3 still 2 — 60 as before. 

3.1 SUGRA19 parameters for low Ahs solutions 

The results from our narrow scan are shown in Fig. [l] as red xs, while points that obey B- 
constraints are labeled as orange xs. The more focused sampling over lucrative parameter ranges 
has produced points with much lower Aew values ranging down to ~ 5, and also solutions with 
Ahs as low as 6. The A^s = 6.4 solution is presented in Tables [T] and |2] as BM model HS2. 
While point HS2 has /i of just 98 GeV, the lightest chargino mass is m^ = 104.1 GeV, slightly 
beyond the limit from LEP2 searches. Since gluino and squark masses are in the several TeV 
range, the point is also safe from LHC8 searches. The mass gaps m^ — m^ = 4.8 GeV and 
m^ — m^ =3.1 GeV so again there will be only tiny visible energy release from the higgsino 
decays. 

To display the sort of parameter choices leading to low Ahs, we show in Fig. |4]the values 
of Aew (blue points) and Ahs (red/orange points) versus superpotential higgsino mass /i from 
the broad (circles) and narrow (x's) scan. From the plot, we see unambiguously that low 
|/i| ~ mz is a necessary, but not sufficient, condition to obtain both low Aew and low Ahs- 
This translates into the solid prediction that four light higgsinos should lie within reach of a 
linear e~^e~ collider with y/s > 2|/i|. 

In Fig. [sj we show Ahs and Aew vs. mH^irricuT) from the broad and narrow scans over 
SUGRA19 parameter space. Here, we see that low Aew solutions can be obtained over a large 
range of mH^iincuT) values, as expected from radiative natural SUSY results[25] which allow 
for a large cancellation between m'jj {rriGUT) and Sm'jj . However, the low Ahs solutions are 
only obtained for rnnS'^GUT) not too far from m^, as required by Eq. [TJ 

In Fig. ^ we show Ahs and Aew vs. M3, where we note that nig ~ IM3I up to radiative 
corrections. Low Aew values allow for M3 ~ 1 — 3 TeV, in accord with LHC searches which 
require rrig ~ 1 TeV for not-too-compressed spectra. 

In the two frames of Fig. [7| we show Ahs and Aew vs. a) M1/M3 and h) M2/M3, where 
all gaugino masses are GUT scale values. The narrow scan has focused on the region around 
M3 ~ 1 TeV, where solutions with Ahs ~ 10 — 100 can be found. From both frames, we see that 
non-universal GUT scale gaugino masses are required for low Ahs solutions, with ratios in the 
range M1/M3 and M2/M3 ~ 2 — 4 being preferred. As remarked earlier, the large electroweak 
gaugino masses provide an initial upwards evolution of m?H which is later cancelled by the 
downward push from the top Yukawa coupling at lower Q values. 

3.2 Sparticle mass spectra from low /Ahs solutions 

We have already seen from Fig. |4] that very low values of /i ~ 100 — 300 GeV are required 
for both low Ahs and low Aew solutions. This translates into the requirement of four light 
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Figure 4: Plot of Ahs and Aew vs. /i from scan over SUGRA19 model parameter space. Color 
coding as in Fig. [ij 
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higgsino states W^ and Z12 with mass ~ 100 — 300 GeV which are difficult to observe at LHC 
but should be observable at a linear e^e~ collider. 

What of gluinos and squarks? We have already seen that rrig ~ 1 — 3 TeV is allowed for 
low Aew solutions while low Ahs solutions prefer a lower value rrig ~ 1 TeV. This is due to 
the fact that too large a value of Ms{mGUT) (and hence a larger physical gluino mass) would 
cause an increase in 3rd generation squark masses, which would increase the Xj factor in the 
m^^ RGE, and cause a larger value of Smj^^ to ensue. 

In Fig. [sl we show A^s and Aew values versus min[mg], where the min is over all physical 
squark masses of the first two generations. While low Aew solutions allow for a broad range 
of rriq ~ 1 — 5 TeV, the low Ahs solutions tend to favor rriq ~ 2 — 3 TeV, beyond current LHC 
reach. 

In Fig. M we show A^s and Aew versus top squark masses a) mi^ and h) mi^. The lowest 
Ahs solutions are found for ini_^ ~ 1 TeV and irti^ ~ 2 TeV. These values are typically 1-2 TeV 
lower than results from RNS models but still beyond most LHC reach projections for third 
generation squark detection. The lowest Ahs solutions, colored dark blue and red, are usually 
in violation of the 5-constraints. This refiects the fact that low Ahs requires light higgsino-like 
charginos (low /i) and light stops, and hence there tend to be large anomalous contributions to 
the BF{b — )■ 57) branching fraction. 

In Fig. [To| we show Ahs and Aew vs. rriA- Here, we see that the lowest A solutions favor 
TTiA (and hence mn and mH±) in the 2-4 TeV range, usually well beyond any projected LHC 
or ILC reach. 
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3.3 Low Ahs solutions and 5-physics constraints 



In Fig. 11, we show values of Aew aiid A^s vs. 



>" 



The recent LHCb 



3.2+1-^ X 10"^ in accord 



BF{B, -^ /i 
measurement [48j finds the branching fraction of BF{Bs — )■ /^^/i^) — 'J-^-i.2 
with the SM prediction of (3.2 ± 0.2) x 10^^. In supersymmetric models, this flavor-changing 
decay occurs through pseudoscalar Higgs A exchange |15], and the contribution to the branching 
fraction from SUSY is proportional to ^ ^^^' . The decay is most constraining at large tan/3 
and at low m^. In the case of low A^s solutions with lower tan/? and heavier m^, we find 
the bulk of solutions to lie within the newly measured error bars although some solutions with 
large Ahs and Aew will be excluded. 

BF{h -^ S7) 



In Fig. 12 



we show Aew and Ahs vs. BF{h — t- 57). We also show the measured central 
value and both 1 and 3-0" error bars. SUSY contributions to the 6 — )■ 57 decay rate come mainly 
from chargino-stop loops and the W-charged Higgs loops, and so are large when these particles 
are light and when tan/3 is large |43t 1^ . In the case shown here, the low Ahs solutions which 
require third generation squarks somewhat heavier than generic Natural SUSY but somewhat 
lighter than radiative Natural SUSY, we find the bulk of low Ahs solutions to lie a couple 
standard deviations below the measured value. 



3.4 Low Ahs solutions and dark matter 

In this section, we show values of the relic neutralino abundance from our scan over SUGRA19 
parameter space. We use the IsaReD[50j relic density calculator from Isajet. Our results are 
frame a), where we plot Ahs and Aew vs. VL^ h"^. For comparison, we also 
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shown in Fig. 

show the location of the PLANCK-measured[CT| relic density of dark matter (green dashed 
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Figure 11: Plot of I^ew and l^us vs. BF{Ba — ^ A*^/^~) from a 19 parameter scan. Color 
coding as in Fig. [TJ The vertical solid line is the measured value and the dashed lines are the 
Icr uncertainties. 
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Figure 12: Plot of I^ew and IS-us vs. BF{b — )• 57) from a 19 parameter scan. Color coding as 
in Fig. [ij The vertical solid line is the measured value and the dashed lines are the la and Scr 
uncertainties. 
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line). Both the low Ahs and Aew solutions populate a band located well below the measured 
abundance. This reflects the fact that the low Ahs,ew solutions all have low fi so that the Zi is 
dominantly higgsino-like; these solutions enjoy an ample annihilation cross section into WW, 
ZZ etc. in the early universe. Thus, the lowest Ahs,ew solutions are typically suppressed by 
factors of 10 — 50 below the measured dark matter abundance. Clearly, additional physics is 
needed in the early universe to gain accord with experiment. One suggestion- the presence 
of late-decaying scalar fields- can either augment or diminish the relic abundance from its 
standard value [52]. Another possibility- mixed neutralino plus axion dark matter- is favored 
by SUSY models with a standard underabundance of neutrahnos since thermal production of 
axinos and thermal/non-thermal production of saxions followed by decays to SUSY particles 
can augment the neutralino abundance [53j (depending on additional Peccei-Quinn parameters 
and the re-heat temperature Tr after inflation). Any remaining gap in the generated neutralino 
abundance relative to the measured abundance can be made up of axions. 

In Fig. [I3^), we plot the rescaled spin-independent neutralino-proton scattering cross section 
{Q~'^h'^/0.11)cr^\Zip) from IsaReS |5l] . The prefactor accounts for the possibility that the 
local dark matter abundance may be well below the value assumed from neutralino-only CDM. 
For reference, we show the latest Xe-100 limit ^55j (dashed green line) for m-^ = 150 GeV. 
The rescaled direct detection (DD) cross section is below current sensitivity by 1-2 orders of 
magnitude. This is lower than even DD projections from radiative natural SUSY[nS] since 
in that case the Zi is mainly higgsino but with a non- negligible gaugino component. Since 
the Zi — Zi — h coupling depends on a product of gaugino-higgsino components, the value 
of a^\Zip) in the RNS case never gets too small. In our current case, with non-universal 
gaugino masses, the Zi can be a much more pure higgsino state, and consequently can have a 
significantly lower DD rate. 

4 Conclusions: 

In previous studies, the radiative natural SUSY model has emerged as a way to reconcile low 
EWFT with lack of SUSY signals at LHC8 and the presence of a light Higgs scalar with mass 
rrih ~ 125 GeV. The RNS model cannot be realized within the restrictive mSUGRA/CMSSM 
framework, but can be realized within the context of NUHM2 models (which depend on 6 input 
parameters) and where /i can be a free input value. In RNS models, Aew as low as ~ 10 can 
be generated while Ahs as low as 10^ can be found. 

In this study, we have implemented scans over the most general minimal flavor- and minimal 
CP-violating GUT scale SUSY model- SUGRA19- with two goals in mind. Our first goal was to 
check if the additional freedom of 13 extra parameters allows for much lower Aew solutions. In 
previous work- by proceeding from mSUGRA to NUHM2 models- a reduction in the minimum 
of Aew of at least a factor of 10 was found[30l ^5\. In the present work, we do not find 
any substantial reduction in the minimal Aew value by proceeding from the NUHM2 model 
to SUGRA19. The parameter freedom of NUHM2 appears sufficient to minimize Aew to its 
lowest values of ~ 5 — 10. 

Our second goal was to check whether the additional parameter freedom can improve on 
the high scale EWFT parameter A^s- In this regard, we find improvements by factors ranging 
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Figure 13: Plot of ^ew and ^hs vs. a) standard neutralino relic abundance VL~h and h) 



rescaled spin-independent o{Z\'p) from a 19 parameter scan. Color coding as in Fig. hi The 
vertical dashed line in a) is the Planck measured value, while in 6) it is the recent Xe-100 upper 
limit on spin-independent neutralino-proton scattering for a 150 GeV WIMP. 
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up to ~ 150! In order to generate low values of A^s, one must generate /x ~ 100 — 300 GeV as 
usual, but also one must start with m|^^ ~ m| at the GUT scale, and then generate relatively 
little change Sm'jj^ during evolution from rncuT to rriweak- Small values of Sm^j^ can be found 
if one begins with electroweak gaugino masses Mi^2 ~ 3M3 at the GUT scale so that gaugino- 
induced RG evolution dominates at high Q ~ rncuT- Then at lower Q values approaching the 
weak scale, top- Yukawa terms dominate the running of m'jj , leading to broken electroweak 
symmetry, but also to not much net change in mjj^ during its evolution from mauT to ruweak- 
The solutions with low A^s are characterized by the presence of four light higgsinos W^ and 
Zi^2 similar to RNS models. However, in contrast to RNS models, the third generation squarks 
tend to be lighter (although not as light as generic natural SUSY which favors m^^ ^ ~ 500 
GeV). The lighter third generation squarks lead to significant SUSY contributions to the decay 
b -^ 57, and seem to be disfavored by the measured value of this branching fraction. In the case 
of low Ahs niodels, the lightest neutralino is more higgsino-like than in RNS models, leading 
to even lower values of predicted relic density and low direct detection rates. The remaining 
CDM abundance may be augmented by scalar field or axino/saxion production and decay in 
the early universe, and in the latter case, the additional presence of axions is expected. 
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